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rn ! Abstract 

^ \ We analyze the phase diagram of single and bi-layer cuprates using the gauge- 

D ' field description of the t-J model. For T > Tbe the in-plane fermion-pairing 

Q ' order parameter Ay is eliminated by gauge field fluctuations, leading us to 

fSJ . predict the absence of a spin-gap phase in single-layer cuprates. For bi-layer 

! cuprates the inter-layer order parameter Aj^ is enhanced by spin correlations, 
and is less affected by gauge fluctuations. We believe that the spin gap in 

' bi-layer materials is due to inter-layer fermion pairing, and that for these 

' underdoped samples the superconducting gap may be nodeless. 

^ ; PACS numbers: 74.20.Mn, 75.10.Jm, 74.25.Nf 
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It is widely accepted that the high-Tc copper-oxide compounds are Mott-Hubbard insu- 
lators which become metallic and superconducting upon hole doping. Thus a useful starting 
point for describing the copper-oxides is the two-dimensional Hubbard model, and many 
workers have studied the strong-coupling limit which leads to the t-J model. A common 
technique to incorporate the constraint of no double occupancy is to decompose the electron 
operator c|g. as the product of a fermion operator f}^ and a boson operator 6,, so that the 
constraint becomes Yla fLfia + bjbi = 1. In this decomposition the fermions carry the spin, 
given by Sj = fl^f^apfip^ and the bosons represent the empty sites. The exchange term can 
be written as 

Si ■ Sj = —\{ficrfja){f]afi<T) ~ + (l) 

= -lU'i-ifji ~ fiifh^U'jifi-i - fnfii), (2) 

which suggests two forms of mean-field decouplings = J^aifLfjir) and Aj^ = 

{filfil ~ fiifi-\)- Let US briefly discuss the mean-field phase diagram, which is quite simple 
for intermediate doping x ^ 0.04 0. It shows the onset of a nonzero below a tem- 
perature T ~ 0.2 J, and the onset of rf-wave fermion pairing (i.e., Aj = — Aj 44.^) at a 
lower temperature shown by the dashed line in Fig. ^ The physical interpretation of this 
state is the formation of a singlet ground state (the RVB state with a gap structure 
in the spin-excitation spectrum. In addition, the bosons are expected to effectively Bose 
condense, i.e. the diamagnetic susceptibility xb will diverge exponentially below a charac- 
teristic temperature, which in the mean-field approximation is given by Tg^ = Airt^x. In 
the region bounded by Tg^ and the onset of fermion pairing, the pairing amplitude {c^Cji) 
of the physical electrons is nonzero, and a (i-wave superconductor is predicted. However, 
above Tg^ the holes are incoherent while the gap structure in the spin-excitation spectrum 
remains, and this region of the phase diagram has been called the spin-gap phase. 

Experimentally there are clear signatures of spin-gap formation below T ~ ISOi^ in 
YBa2Cu306+x with x = 0.6 and in YBa2Cu408. The signature is clearest in the copper and 
oxygen nuclear spin-relaxation rates ^^Ti^ and ^^T]~^, but it is also apparent in the Knight- 
shift data and in T2 measurements 0. Since the spin gap appears only in underdoped 
samples, one is naturally led to an explanation in terms of the fermion pairing discussed 
above. However, recently Millis and Monien re-examined the experimental evidence , and 
concluded that the spin-gap phase is absent in single-layer material such as La2_xSra.Cu04. 
They suggested that the inter-layer exchange couphng J± may be responsible for the spin- 
gap phase. However in the models they constructed, an unreasonably large exchange J±, 
comparable to J, is required to explain the experimental data. 

On the theoretical side, there has been much effort on improving the mean-field theory 
by including fluctuation corrections. It has been recognized that the phase fluctuations of 
the field are the most important, leading to a gauge theory in which the fermions and the 
bosons are coupled to a U(l) gauge field PJ^. Integration over the gauge field enforces 
the condition that a flow of fermions is accompanied by a backflow of bosons. Initially the 
gauge field has no dynamics, which it only acquires upon integrating out the fermions and 
bosons. We focus on the transverse gauge field because, unlike longitudinal fluctuations, 
transverse currents are not screened, giving rise to unusual infra-red singularities 0. The 
transverse gauge-field propagator is given by D(q, cu) = H(q, a;)~^, where H = H^ -|- H^, 
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ni^(q, iun) = XF<f + It^nl^F/*?, and nB(q, iuJn) = XB<f - Here Xf and Xb are the fermion and 
boson diamagnetic susceptibility, and for free bosons Xs = 

(exp(T]^°i^/T) - l)/247rmij, while 
the \uJn\kF/q term represents Landau damping. Recently we found that gauge fluctuations 
modify the mean-field phase diagram in a significant way P,^. The Bose-condensation 
temperature is greatly suppressed, and the spin-gap phase is eliminated entirely. The former 
effect is due to thermal fluctuations of the gauge field, given by (h(x) ■ h(0)) ~ (T/x)(5(x) 
where h = V x a and X = Xf + Xb- Note that at sufficiently low temperatures the RMS flux 
per plaquette is small. Nevertheless, since the deBroglie wavelength At of the bosons covers 
many lattice spacings, the effect on the random gauge flux is strong. Indeed we find that the 
dimensionless coupling constant \s g = (T/x)^t ~ {™bx)~^- At temperatures above Bose 
condensation we may ignore x'b compared to Xf that g = 247rt/ J, i.e. the system 

is in the strong-coupling limit even for modest t/J. In this limit the Feynman paths of the 
bosons are almost self-retracing. By taking only static gauge-field configurations and taking 
a quenched average, we find that xb grows inversely with T, but is suppressed compared to 
the high-temperature limit of Xb ^ factor {0.2g)~^ 0. At sufficiently low temperatures, 
Xb becomes much larger than x^f ^^'^ "^^^ therefore no longer be ignored in the estimate of 
g. We crossover to the weak-coupling limit and xb self-consistently diverges exponentially 
below a temperature Tbe which we estimate to be This suppression of the effective 

Bose condensation eliminates one of the most serious drawbacks of the mean-field phase 
diagram, where T^^ yields a temperature scale which is much too large. 

We also examined the effect of gauge fluctuations on the fermion pairing. Here we find 
that it is the quantum fluctuations which are important, and the effect is strong enough to 
de-stabilize the mean-field transition until the temperature falls below Tbe [H- We computed 
the gauge field contribution to the free energy 

^ sr^ f°° du , , ^, flmIl(q,LJ + iS)\ 

"^^-- = ^70 ^[2^BM + l]arctan(^ ^^^^^^^^^^J . (3) 

For T > Tbe, H may be replaced hj Up, and Eq. (|^) yields a large negative contribution to 
the free energy which is cutoff by T, leading to a T^^'^ correction term, and a specific heat 
of T^/'^ . Now we can ask what happens if the fermions form a d-wave pairing state with 
amplitude A. We can compute TgaugelA) by replacing H in Eq. (^) by the corresponding 
polarization function H(q, ia;„. A) in the pairing state. Since a gap appears in ImH(q, cu), 
the u integral is cutoff by A instead of T, and we obtain Fgauge ~ A^/^. This is to be 
added to the mean- field gain in free energy Tmf ~ — A^, which is overwhelmed by Tgauge for 
small A. However, a first-order transition is in principle still possible. We have carried out a 
detailed numerical computation of Tgauge(A) 0, and found that the pairing state is unstable 
for T > Tbe- Below Tbe, Xb{T) grows exponentially and consequently Re H^ dominates 
over Re Hi?. Thus FgaugelA) becomes less important and the cost of opening up a gap is 
correspondingly weaker. Our calculations show that a first-order transition to a ci-wave 
pairing state occurs below Tbe- Since the bosons are coherent, this is a direct transition 
to a d-waye superconductor. Our best estimate of the transition temperature is shown in 
Fig. |1|. We expect the first-order nature of the transition to be smoothened out once phase 
fluctuations of the superconducting order parameter are taken into account. An important 
consequence of this analysis is that the spin-gap phase is completely eliminated by gauge 
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fluctuations. This agrees with the absence of a spin-gap phase in single-layer materials such 
as La2_xSrj.Cu04. 

Now that we have de-stabilized the spin-gap phase in a single layer, we have to appeal 
to inter-layer coupling to explain the appearance of spin gaps [Q . In a bi- layer material it 



is natural to look for an order parameter 10 



Ax(r,) = {/,';'/if-/i"/if), (4) 

by decoupling the inter-layer exchange term = J±J2iSi^^Sf^ . At present we only have 
a lower bound J± ^ 8 meV from the absence of an optical mode in the neutron scattering 
from undoped YBCO |ll]], but we expect Jj_ to be much smaller than J. A conventional 
pairing theory would then predict a gap which is exponent ionally small in J/Jj_. To over- 
come this we note that the underdoped cuprates have a reasonable long antiferromagnetic 
correlation length in the plane. In a bi-layer it is then favorable to form singlets not only 
between spins which are directly on top of each other, but also between other pairs on the 
same sublattice. Thus we expect the inter-layer pairing gap to be enhanced by in-plane 
correlations. To analyze this we employ the random-phase approximation (RPA), which has 



been quite succesful in reproducing the experimentally observed spin correlations 0]. The 
RPA approach corresponds to a decoupling of the exchange term Sj ■ S-, in a third channel 
which differs from the decouplings in Eqs. (P and (|^). In order to avoid double counting 
we separate the Hamiltonian into two equal parts, decouple the first part using Eq. (|1|) and 
apply RPA on the second part. Our results should not be too sensitive to this admittedly 
ad hoc procedure. Within this scheme we define Jf^ici) as the effective inter-layer exchange 
due to the summation of a series of RPA bubbles shown in Fig. 1^. This analysis will be 



discussed in more detail in Ref. [1131. The result is 



where x'^(q) is the free-fermion bubble for the two-dimensional band structure. The de- 
nominator vanishes at the onset of the Neel ordering, which occurs at Xc — 0.08 within our 
approximation. The numerical solution of Jj^{c{) for shown in Fig. ^ and we see 

that it is significantly enhanced near the incommensurate nesting vector Qaf — (vr, vr ± 2x). 
We find that for x = 0.09 the peak at Qaf leads to an antiferromagnetic correlation length 
of approximately 3 lattice spacings. We next take the Fourier transform of Jx^(q) to obtain 
Jf^{rij) and consider the effective inter-layer exchange J^^(rjj)S,-^'' ■ S^-^"*, which is decoupled 
using Eqs. (||) and (|). The gap equation for the order parameters A±{rij) is solved numer- 
ically in real space. We find that the pairing amplitude is especially large on the diagonal 
such that Aj_(rij) ~ |Aj|(— 1)* decays slowly with i. Transforming A(rjj) back to k space, we 
find a quasi-particle dispersion E(k) = [e(k)^ + Ax(k)]^/^, where A_L(k) is shown in Fig. |^. 
Notice that Aj_(k) exhibits some anisotropy but does not change sign. Thus the inter-layer 
pairing is s-wave in nature, with a full gap in the excitation spectrum. 

It is important to note that it is crucial in this approach that the in-plane rf-wave pairing 
is suppressed by gauge fluctuations. One might ask whether gauge fluctuations will destroy 
the inter-plane pairing as well. We believe that this does not happen for the following 
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reason. There are two gauge fields modes ai and 02 in the two layers, which are split into 
a± = (ai ± a2)l\/2 at the onset of inter-layer pairing. The gauge propagators are given by 



n_|_^ where 



n±(q,zi/0 = c + 2r5] 1 



i2n) 



ee' - UnUj'n ± AA' 




(6) 



where iu'^ = iun - iz/„; e, e' = e(k ± q/2); A, A' = Ax(k ± q/2); and E = Ve^ + A^. For 
A = 0, the constant C exactly cancels the second term for q ^ and Un = 0, so that 
the gauge field is massless in the normal state. In the pairing state a gap opens up in 
n+(o,o) ~ Ai, but n_ remains massless. Furthermore, when we substitute Il± into Eq. 
(0), we find that while the 11+ propagator introduces an energy cost and is pair-breaking as 
before, the n_ propagator is actually pair- enhancing. This is because even though an energy 
gap is introduced into Im n_, the coherence factor in Eq. (|^) is such that — (Imn_/Ren_) 
is enhanced for u > 2A± compared to its normal state value, and this overwhelms the loss 
of the contribution from uj < 2A±. This pair-enhancing nature of the a_ mode can also be 



understood in another way |jTj] . The fermions on the two planes couple to the a_ mode with 
opposite charge, so that the exchange of an a_ mode leads to an attraction, analogous to 
what happens in the t-t'-J model [0. In our case we expect that the effects of the a+ and 



a_ gauge fluctuations largely cancel each other, leaving the mean-field transition intact. 

Since we are unable to compute A^ quantitatively, in the remainder of the paper we 
will use A^ as a parameter to compute various physical quantities, which we compare with 
experiments. We calculated the nuclear- relaxation rate 1/Ti for the copper and the oxygen 
sites, and we immediately encountered the problem that for s-wave pairing the Hebel-Slichter 
peak appears below Tc. We circumvented this by arguing that since the order parameter is 
not invariant under the local gauge transformation /j^'^'* — > e*^'/fj'^'' the transition should be 
broadened into a crossover behavior, and a BCS-type gap should be replaced by a pseudo gap 
as shown in the inset in Fig. |^. With this and a modest broadening F = O.IT, we eliminated 
much of the Hebel-Slichter peak, and the results are shown in Fig. Notice that above the 
crossover transition, (^^TiT)~^ increases upon lowering T, which is due to antiferromagnetic 
correlations. We do not see this increase in (^^TiT)~^, because the form factor ^'^F(q) for 
the oxygen site vanishes at the antiferromagnetic nesting vector Qaf • We also computed 
I/T2 and the uniform susceptibility x(q = 0), shown in the second panel of Fig. ^ The 
results compare reasonably well with experiments provided we choose Aj^ — 150K ~ O.IJ. 
We believe that this is not an unreasonable value if J± ~ 0.2J. 

So far we have ignored any inter-layer hopping of the form t±c\^'^c^-^ , which is reasonable 
provided xt^ < J±. If this is violated we expect inter-layer pairing to be suppressed, but we 
have not studied this quantitatively. Not enough is known about t± and J±, but our guess 
is that xt± and J± are comparable. However, even a small t± will lead to coherence between 
bosons on the two planes immediately below Tbe, so that the fermion pairing becomes 
genuine superconducting pairing between electrons on the two layers. At low temperatures 
the in-plane s-wave and the inter-plane d-wave pairing will co-exist, giving rise to a quasi- 
particle dispersion E{k) = {e(kf + A±(kfy/\ where A±(k) = A^(k) ± A||(k). If A^ 



5 



is indeed as large as 150 K, as the experiments seem to indicate, comparison with Fig. |l| 
indicates that for underdoped materials it is likely that A_|_ > |A||| for all k, so that the 
superconducting gap is nodeless. As doping is increased, Aj^ decreases rapidly with the loss 
of antiferro magnetic correlations, which was essential for the enhancement of J^, and we 
crossover to a superconducting state with nodes. As long as A_l remains finite we expect that 
there are eight nodes along the Fermi surface instead of the four nodes for a conventional 
d-wave superconductor. This new phase diagram is shown schematically in the insert in Fig. 
|l|. Finally we note that Chakravarty et al. have shown that inter-layer pair tunneling 



leads to an anisotropic s-wave gap function. We emphasize that in our theory the order- 
parameter symmetry is always d-wave in the plane, even in the nodeless state. Furthermore, 
to the extent that our nodeless pairing state is always preceded by a spin-gap phase, we 
believe that the origin of superconductivity in our treatment is closer to the original RVB 
picture of pre-formed spin-singlet pairs [|l| than the more recent inter-layer pair-tunneling 
model. 

We acknowledge helpful conversations with A.J. Millis. This work was supported by the 
NSF through the Material Research Laboratory under Grant No. DMR-90-22933. 



6 



REFERENCES 

[1] P.W. Anderson, Science 235, 1196 (1987); G. Baskaran, Z. Zou, P.W. Anderson, Solid 

State Commun. 63, 973 (1987). 
[2] G. Kotliar, J. Liu, Phys. Rev. B 38, 5142 (1988); M.U. Ubbens and P.A. Lee, Phys. 

Rev. B 46, 8434 (1992). 
[3] See H. Monien, D. Pines, and M. Takigawa, Phys. Rev. B 43, 258 (1991). 
[4] A.J. Millis and H. Monien, Phys. Rev. Lett. 70, 2810 (1993); E 71, 210 (1993). 
[5] G. Baskaran and P.W. Anderson, Phys. Rev. B 37, 580 (1988); G. Baskaran, Phys. Scr. 

T 27, 53 (1989). 
[6] L.B. loffe and A.I. Larkin, Phys. Rev. B 39, 8988 (1989). 
[7] P.A. Lee and N. Nagaosa, Phys. Rev. B 46, 5621 (1992). 
[8] M.U. Ubbens and P.A. Lee, M.I.T. preprint (1993). 
[9] M.U. Ubbens and P.A. Lee, to appear in Phys. Rev. B (1994). 
[10] B.L. Altshuler and L.B. loffe. Solid State Commun. 82, 253 (1992). 
[11] J.M. Tranquada et al, Phys. Rev. B 46, 5561 (1992). 

[12] T. Tanamoto, H. Kohno, and H. Fukuyama, J. Phys. Soc. Jpn. 61, 1886 (1992); J. Phys. 

Soc. Jpn. 62, 717 (1993). 
[13] M.U. Ubbens and P.A. Lee, M.LT. preprint (1993). 
[14] We are grateful to N. Bonesteel for a discussion on this point. 

[15] P.B. Wiegmann, Phys. Rev. Lett. 60, 821 (1988); Physica (Amsterdam) 153-155C, 
103 (1988); X.G. Wen, Phys. Rev. B 39, 7223 (1989); P.A. Lee, Phys. Rev. Lett. 63, 
680 (1989). 

[16] S. Chakravarty et al, Science 261, 337 (1993). 



7 



FIGURES 



FIG. 1. The phase diagram for single- layer cuprates. The dashed line is the mean-field pairing 
transition. The solid line takes gauge fluctuations into account, so that fermion pairing only survives 
for T < Tbe, resulting in a d-wave superconductor. The inset is a schematic phase diagram for 
bi-layer materials. The spin-gap phase is due to inter-layer pairing. The superconducting state has 
d symmetry in the plane with nodes in the gap function, but there may be a region indicated by 
the shaded area where the gap is nodeless. 

FIG. 2. The effective interlayer coupling J^{ci) for various values of the doping x, using 
J± = 0.2 J. This is obtained by summing over RPA bubbles (inset). For x — 0.08 there are strong 
incommensurate peaks at the nesting vector Qaf- 

FIG. 3. The gap Ax(k) for x = 0.085 and Jj_ = 0.2 J. The gap has an extended s-wave 
symmetry, and is enhanced close to the Fermi surface (dotted line), especially at the corners. 

FIG. 4. The first panel shows the nuclear-relaxation rates (TiT)~^ on the ^^Cu and the ^^O 
sites. The rise is (^^TiT)^^ is due to antiferromagnetic correlations, which is absent in (^^TiT)"^. 
The second panel shows I/T2 and x(0, 0). The opening of a pseudo-gap (see inset) causes a rapid 
decrease in (TiT)''^, I/T2, and x(0,0) for T < O.IJ. 
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